Nonlinear bending of molecular films by polarized light 
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A theory of photoinduced directed bending of non- crystalline molecular films is presented. Our 
approach is based on elastic deformation of the film due to interaction between molecules ordered 
through polarized light irradiation. The shape of illuminated film is obtained in the frame of the 
nonhnear elasticity theory. It is shown that the shape and the curvature of the film depend on 
the polarization and intensity of the light. The curvature of an irradiated film is a non-monotonic 
function of the extinction coefficient. 

PACS numbers: 46.25.Hf, 46.70.De, 68.55.-a, 78.20.Bh 



I. INTRODUCTION 

Polymer films and solids containing light-sensitive 
molecules have the remarkable property to change their 
shape and size when irradiated with light. Certain poly- 
mer films containing azobenzene chromophores in the 
main and side chains exhibit strong surface relief fea- 
tures under illumination: trenches under the action of 
linear polarized light and mounds or wells under the ac- 
tion of circularly polarized light [11, [3] . Circular azoben- 
zene polyester films freely lying on a water surface be- 
come elliptically deformed under the influence of linearly 
polarized light [3'] . Large, reversible shape changes can be 
induced optically by photoisomerization of nematic elas- 
tomers [4] . Anisotropic bending and unbending behavior 
of molecular liquid-crystalline films containing azoben- 
zene chromophores has been discovered and studied in 
Refs. d, [1] where it was shown that the films can be re- 
peatedly and precisely bent along any chosen direction 
by using linearly polarized light. Fast (on the timescale 
of 10~^ s) light induced bending of monodomain liq- 
uid crystal elastomers has been observed observed in 
0. Shap e-memory effects in polymers containing cin- 
namic groups induced by ultraviolet light illumination 
were reported quite recently in Ref. Q. The possibility 
of coupling between orientational and translational de- 
grees of freedom in liquid crystals was first raised b y d e 
Gennes @ and extended to nematic elastomers in [lC| . 
Based on this idea a phenomenological theory of pho- 
toinduced deformations of nematic elastomers was pro- 
posed in [ll| . A microscopic theory of photoinduced de- 
formation of non-crystalline molecular films was devel- 
oped in [T^ . The physical reason for surface-relief for- 
mation presented in is that azo-dyes have two iso- 
meric states: cis and trans. The molecules in these two 
states have significantly different shapes. For example, 
in the case of azobenzene chromophores the irans-isomer 
is highly anisotropic whereas the cis-isomer is approxi- 
mately isotropic [l3[ , so the multipole moments and sizes 
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may differ significantly. It was shown that there are 
two contributions to the photoelastic interaction: from 
the orientational interaction between molecules and from 
the interaction which is due to the change of the van 
der Waals interaction energy between a molecule and all 
surrounding molecules in its transition to the cis-isomer 
state. The former causes the film deformation under 
the action of linearly polarized light while the latter (to- 
gether with the orientational interaction) is responsible 
for the surface relief formation under the action of circu- 
larly polarized light. The possibility of creation of wells 
and humps on the film surfaces under the action of cir- 
cularly polarized light was discussed. It is worth not- 
ing that in the frame of this approach neither orienta- 
tional order nor orientational in-plane anisotropy of the 
film in the absence of irradiation was assumed. One can 
say that the absorption of linearly polarized light cre- 
ates an orientational order in the film which in turn pro- 
duces anisotropic deformation. Quite recently, based on 
the idea of isomeric states, a theory of the polarization- 
dependent photocontractions of polydomain elastomers 
due to light-induced director rotation, was proposed in 

The aim of the present report is to apply this approach 
to anisotropic bending of molecular films by polarized 
light. We present an elastic energy of the film in the 
presence of polarized light and solving equations of equi- 
librium, we show that a change in the polarization direc- 
tion of light causes a corresponding change of the shape 
of the film. We also show that the curvature of an irra- 
diated film is a non-monotonic function of the extinction 
coefficient. 



II. ELASTIC ENERGY OF MOLECULAR FILM 

We consider a film containing molecules with two dif- 
ferent isomeric states. Let the middle surface of the film 
coincide with the x,y-plane so that the undcformcd film 
occupies the region: \z\ < x,y ^ The film is irradi- 
ated by a linearly polarized electromagnetic wave which 
propagates along the z-axis. Its electric component has 
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the form 



E{r, t) = £ cos {k z — ujt) 



(1) 



where lo is the frequency, k — to/c is the wave number (c 
is the speed of hght) and 



£ = £ (cos sin ip, 0) 



(2) 



is the amplitude. The angle '0 determines the polariza- 
tion of the electromagnetic wave. 

The total elastic energy of irradiated thin film may be 
written as follows 



(3) 



Here the first term represent the elastic energy of the 
non-irradiated film 
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where E and a are Young's modulus and Poisson's ratio, 
respectively [l^ and ea p is the strain tensor (a, (3 — x,y). 
We assumed here that that the non-irradiated film is nei- 
ther orientationally nor translationally ordered in the xy- 
plane and therefore we modeled its elastic properties by 
using the isotropic energy The term 



Wg = P I dxdy w 



(5) 



in Eq.([3]) presents the potential energy of the film in the 
gravitational field. Here P = agr pf his the gravity force 
with p f being the film density and agr being the accelera- 
tion of free fall. The last term in Eq. ([3]) gives the change 
of the elastic energy due to the interaction between the 
electromagnetic wave ([T|) and the film. It has the form 

m 



W = Wi + W2 , 



(6) 



where 
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Wi = -Va J J dxdy N{r) {e^x - eyy) cos{2tp) 

-h/2 O 

+2exySin{2iP) (7) 

represents the anisotropic part of the photoelastic inter- 
action which describes the coupling of the shear defor- 
mation of the film to the incoming light and 

h/2 

W2 = -Vi j dzJdxdyM{r){exx + eyy) (8) 
-h/2 n 



determines the isotropic in-plane deformations. In 
Eqs.([7])-(I5]) the parameter of the photoelastic interaction 
Va is due to the orientational (e.g., dipole-dipole) part of 
intermolecular interaction while the parameter Vi is due 
to isotropic part of the intermolecular interaction. For 
the sake of simplicity we assume them space independent. 
The function J\f(f) gives the population of cis-isomers 
for a given value of the radiation power £^ J2]. Fol- 
lowing the Bouguer-Lambert-Beer law which determines 
how the intensity of light decreases under its propagation 
inside an absorbing medium, we shall model the function 
7V(r) as follows 



Mir) = Ma cxp 
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where A/o is the maximum population of the czs-isomers 
for a given power £'^ and ^ is the extinction length of 
the light which provides transition of chromophores from 
trans- to cis-isomeric state (in the experiments it 
was the light with the wavelength 366 nm); other distri- 
bution can of course be used, depending on the actual 
arrangement of azo-dyes. Note that in Eq.® we ne- 
glected the fact that upon bending the normal to the film 
surface deviates from the z-direction which is legitimate 
when bending is small. 

Note that in the case of circularly polarized light when 
instead of Eqs.dl]), ^ we have 

E{f, t) ^ £ {cos(fcz - ujt) , sin{kz - cut) , 0} (10) 

and therefore the contribution ([7]) vanishes; the energy of 
photoelastic interaction is solely determined by Eq.®. 

Following the usual derivation of Foppl-von Karman 
equations for bending of a thin plate (see, e.g., |il6i]). 
we write the strains as a linear expansion in z from the 
middle plane and get 
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^al3 = ^ {dfSUa + daUjs) + ^ daW dfjW 



(11) 



(12) 



are the components of the two-dimensional nonlinear de- 
formation tensor and da denotes differentiation with re- 
spect to the coordinate x^ = y. Referred to these 
coordinates, the components of displacement are Uq — 
Ux,Uy,w where we have named the vertical displacement 
Uz = w. Introducing this equations into Eq.Q we can 
represent the elastic energy of an irradiated film in the 
form 



F ^ Fb + Fs + W . 
The bending energy is written as 

D 



(13) 



Fb = 



dxdy^ {AwY 
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dlwdlw 
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where A = da da is the two-dimensional Laplace operator 
and 



The function 



D = 



Eh^ 



12(l-cr2) 



(15) 



is the flexural rigidity of the film. For the stretching 
energy one has 



-Fs = 2 y ^2;dy {a^^ u^x + cfyy Uyy + 2a^y Uxy) (16) 
n 



with the longitudinal stresses 
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~ 1 + a ' 
In Ea. p3p the light-film interaction energy becomes 

W = Wb + Ws (18) 
with the bending contribution 

Wb = Wbi + Wb2 , (19) 

where 



Wbi 



dxdy^ 



AiAw 



{dlw - dlw) cos(2V') -t- 2 d^yW sin(2i/') | (20) 



describes a linear interaction with bending deformation 
while the term 



Wb2 = - 2 y dxdy[Ba [( {dxwf - [dywf ) cos(2i^) 
+2 sin{2ip) -f 



;(9,H' + (9,«;)']} (21) 



is due to the nonlinear character of the two-dimensional 
deformation tensor (|12p. The stretching contribution in 
Eq.dlll) is 

= -h j dxdy^Ba {u^x - Uyy) cos(2V') 
n 

-1-2 Ua;^ sin(2V') + Bi{uxx + Uyy)^ , (22) 

and the new light-film interaction parameters in 
Eqs.(l20l)-(I22]) are 

Ba=VaJl. (23) 
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gives the mean value of cis-isomers in the film and 
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(25) 



characterizes asymmetry of the cis-isomer distribution in 
the film. 

By using Green's formula for the two-dimensional inte- 
grals, Eas. l|20| and ((22l) can be presented in an equivalent 
form 



Wbi = -h^ j>dl\^ [A.,+ AaCos2{^ - 9)] 



dw 
dn 



+AaSm2{i;-9)^} , (26) 



W. = -hfdl{B.^n.u + B^cosi2^-0)ux 

+Basm{2i:-e)uy'j , (27) 

where d/dl is the derivative along the tangent I to the 
contour and it has together with the normal derivative 
d/dn the form 
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(28) 



where the angle 9 determines the direction of the out- 
ward normal to the boundary contour: n = (cos 0, sin 0). 
Eqs.([20| and ((22|) present two physically different mech- 
anisms of film deflection. Eq. describes the light- film 
interaction which causes the change of the film area. The 
intensity of the interaction is proportional to the mean 
value of cis-isomers in the film J\f and as it is seen from 
Ea. (P7|) . the action of light is equivalent to a uniformly 
distributed edge force applied in the plane of the film. In 
the presence of the interaction with Bi > the area 
of the film increases while the opposite sign corresponds 
to the compression of the film. On the other hand the 
interaction given by Eq. (j20p is due to asymmetric distri- 
bution of cis-isomers in the film As Eg. ip^ shows, 
in this case the light produces a bending moment applied 
to the boundary contour. 
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III. DEFLECTIONS UNDER THE ACTION OF 
POLARIZED LIGHT 



The Euler-Lagrange equations for the functional ([T3|) . 
(dH), (HH), ([TH])-(I22]) (Foppl-von-Karman equations) have 
the form 



d 



dw 



D w — h — — aaB i- hBiAw 

ox 13 OXa 

+hBa [{dlw - d^w) COS 2-0 + 2 dxyW sin 2ip\ 

= -P , (29) 

-^a^p^O. (30) 

OXfj 

The boundary conditions for these equations may be ob- 
tained in the same way as it was done in [Tsj and in the 
case of free boundary (the edge of the film is free) the 
variations of the vertical component 5w and its normal 
derivative d{6w) jdn on the edge are arbitrary. This gives 
the following set of equations 
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cos20dxyw} - /iMa^sin2(V' - 6*) = , (31) 
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Aw -{I- cr){sin2 edlw + cos^ Bd^w 



-sm29da;yw} ~ h'^[A, + Aa cos 2{ip ~ 9)] = , (32) 

o'xff — {Bi + Ba cos 2tp) cos 9 
—Ba sin2'!/; sm9 — , 

f^vpnp — [Bi — Ba cos2V') sm9 

-Ba sm2iP cos9 = . (33) 

By introducing the Airy potential x(a;,y), so that 
Eas. (|5D)) are automatically satisfied: 

<^xx ^dyX , f^yy^dlx, CTxy = CTyx = "O^y X , (34) 

and presenting the Airy potential as a sum 

X{X, y) = ^ (Bi ~ Ba COS 2^) 



+ -{Bi + Ba cos2V')y^ 
'Ba sm2ipxy + T{x,y) , 



(35) 



we obtain the Foppl-von-Karman equation (|29p in the 
form 

DA'^w - h {dlTdlw + d^T d^w - 2dxyT d^yw) 

= ~P . (36) 

This equation is to be completed by the compatibility 
condition llal 



A^T+ E 



dlwdlw-{dxywY 



= . 



(37) 



Note that owing to the fact that we extracted from the 
Airy potential x(a;,y) the parabolic contribution, the 
Foppl-von-Karman equation p6p does not contain light- 
induced driving terms. Introducing Eqs. ([M)l . ([55]) into 
Eqs. (|33p . we obtain that the boundary conditions for the 
potential T{x, y) have particularly simple form: 



a^T sin6l - d^yT cos9 = , 
dl T cos 61 - d^y T sin 6* = . 



IV. BENDING OF CIRCULAR FILM 



(38) 



We consider the photoinduced bending of a circular 
film of radius R. In the polar coordinates 



x = pR cos 9 , y — pR sin ( 



(39) 



where p — r/R, [r = \J x^ + y^) is a dimensionless radial 
coordinate and 9 is the azimuthal angle we obtain from 
Eas. (|55)) and ([57)1 that equations for new dimensionless 
dependent variables 



w h 



(40) 



take the following form 

A'C - ^{ (a|C + pdpC) + {dlr + pd,T) die 

[deT - pdper) {dgC, - pdpeQ | 

= -/ , (41) 

12(l-a2)-l{(a,2c + p9pC) d^C 

-4 (9eC-p9,eC)'}-0. (42) 



The boundary conditions PT|) . ([5^ and 
(r — R) become 



at p = 1 



dpAC 



-aa cos2(V' -9)^0 



- [a, + aa cos2('(/; - 9)] = 



(1 - pdp) dgr^O, 
pdpT + dlr^Q, 



(43) 



(44) 



(45) 



where / = PR'^ /{hD) is a dimensionless gravity force 
and dimensionless parameters 
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hR^ 



(46) 
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FIG. 1: The shape of the circular film illuminated by a circu- 
larly polarized light. 



characterize the intensity of the light-film interaction. 

To have some insight we consider first the case of weak 
light-film interaction: < 1, Oa < 1. Then we compare 
these results with numerics. 
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A. Circularly polarized light 

In this subsection we consider the case of photoin- 
duced film deformation when the light is circularly po- 
larized pO| (the parameter Oa is set to zero). Assum- 
ing the radial symmetry of solutions [dgC, = der — 0) 
from Eqs. (|4T]l - (|45|) one can obtain approximately (see 
Appendix for details) that the shape of an irradiated film 
is determined by the expression 

(see Fig.ll]). 

In order to verify the range of validity of the ap- 
proximate solution (|47p . full numerical simulations of 
Eqs. ([lT|) - (|45p have been performed for the case of ra- 
dial symmetry. We used finite difference method solving 
the resulting set of nonlinear algebraic equations by New- 
ton iterations. The results of numerical calculations for 
the shape of irradiated film ^(p) are shown in Fig. [5] for 
the different values of the dimensionless gravity force / 
and the parameter ai together with the approximate de- 
pendence (|47p . Even for large value of parameter ai the 
relative difference between numerical solution and the ap- 
proximation (|47)) does not exceed 25% for ai — 10 and 
becomes smaller for smaller Oj. For ai — \ the numeri- 
cal and approximate solutions are undistinguished in the 
Fig.H 



FIG. 2: Dimensionless displacement C,(p) of the circular film 
illuminated by a circularly polarized light. Numerical (solid 
lines) and approximate (dashed lines) solutions for a — 0.5 
and / = 1 (a), / = 10 (b), and different values of a^. 

B. Linearly polarized light 

Let us consider now the case of photoinduced film de- 
formation caused by the linear polarized light ([T]), ([2]). 
Expanding solutions of Eas. (|lT|) - (H51) in a Fourier series 

oo 

C(p,0) =^ Cn(p) cos[2n{^-9)] , 

oc 

r (p, 0) = ^ r„ {p) cos [2 n ( V - 0)] , (48) 

and assuming the weak light-film interaction < 1, Oq < 
1 [the condition Oq < 1 allows to take into account only 
the zeroth and the first components in the Fourier series 
dH])] one finds 

dp, 0) ^ Co(p) - aa ^e(3^y ^ cos2(V^ - 0) (49) 

for the vertical component of the displacement (see 
Fig. [3]). 

For comparison full numerical simulations of Eqs. (j4ip - 
(j45p have been performed by use of finite difference 
method. The results of calculations are shown in Fig. U) 
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FIG. 3: The shape of the circular film illuminated by a linearly 
polarized light. 



For small values of parameter the high order terms 
Cn, Tn in the expansion decrease for n > 1 and the 
numerical solution for ^(p, 6) practically coincide with 
approximation (|49p . 



V. BENDING OF RECTANGULAR FILMS 

Let us now consider the case of rectangular film: x G 
{—L/2, L/2), y G (— L/2, L/2) interacting with a linearly 
polarized light. We will consider the case of weak light- 
film interaction and use the linear theory (neglect the 
stretching contribution). We will also neglect the action 
of the gravity force. As a solution to Eq. ([29|) we choose 



1 2 
w = -ax 



bxy + ]^cy^ 



(50) 



where a, h, and c are some constants. These constants 
can be found by using the boundary condition ([5^ . How- 
ever a more simple way to find them is to introduce 
Eq.dSOl) into Eqs.dHI) and (HH) and get 

l.F=^[{a + cf + 2{l-a){b'-ac)] 
-/iMj(a + c) -/i^a [(a- c) cos 2^^-1-26 sin 2-0] . (51) 
The function (15111 has a minimum for 



a = h 
c = h 



2 A, (l-g) + (l + g) Ag cos2-i/; 
i?(l~a2) 

2 Aj(l - Cr) - (1 + (t) Aa C0S2V' 



Dil 



b = h^ 



Aa sin 2-0 



(52) 



D{l-a) 

Thus under the action of linearly polarized light H]) , ^ 
an initially flat molecular film takes the shape 

w — — Ki {x cos ^ + y sin ■0)^ 
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FIG. 4: Dimensionless displacement (^(p) of the circular film 
illuminated by a linearly polarized light for ?/> = in the cross 
sections — and 9 = 7r/2. Numerical (solid lines) and 
approximate (dashed lines) solutions for a = 0.5 and / = 10: 
= ffla = 1 (a), tti = Oa = 5 (b), and ai = aa = 10 (c). 



+ 2 '*2 {—X sin-0 -I- ycosip)"^ 



(53) 



which is characterized by the following two principal cur- 
vatures 



Ki = h 
K2 = h 



2 A,{1 - a) + Aail + cr) 



2 A^jl - a) - Aa{l + a) 



(54) 
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FIG. 5: The shape of the rectangular film with Aa < Ai 
illuminated by a linearly polarized light for t/) = 0. 

and has the equilibrium energy given by the expression 

F - -^^0^) [^'(1 -^) + AlH + t)] . (55) 

To characterize the global shape of the film it is conve- 
nient to introduce the mean curvature H = [ni + K2)/2 
and the Gaussian curvature K = ki K2. 

When K > Q the point {x — Q,y — 0) is an elliptic one 
and the film has a paraboloid shape. For if = it be- 
comes cylinder-like. Let us assume first that > |k2|, 
then an irradiated film takes a shape close a cylindric 
one. Note that this is probably the case in the exper- 
iments of Refs. @] where figures show that polymer 
films are bent in a cylinder-like fashion. In this case 
when the light is polarized along the x-axis (V' = 0) the 
bending occurs around the y-axis (see Fig.[5l). Under the 
action of light polarized along one of the film diagonals 
(-0 = 7r/4, 37r/4) the film bends around a diagonal (see 
Fig. [HI). This behavior is in a full agreement with the 
results of the Refs. @]- When K < Q the film takes a 
saddle-like shape and the corresponding shape profile is 
shown in Fig. ([7]) . Note that this probably the case in the 
experiments of Ref. Q where such of kind of deformation 
was observed for liquid-crystal elastomers with azo-dyes. 

The principal curvatures of the film ki and k,2 are lin- 
early proportional to the maximum population of the 
cis-isomers in the film A/q and they are a non-monotonic 
function of the extinction length In Fig. [5] we presented 
the normalized mean curvature H / Hm {Hm is the max- 
imum value of the mean curvature) as a function of ^, 
restricted to ^ < 2/i {h is the film thickness). The last 
statement may be considered as a microscopical explana- 
tion of the conclusion which was drawn in Ref. |JJJ in the 
frame of a phenomcnological theory. Note also that it is 
clear from Eqs. ((23)) . ((25)l and (|54)l that for a decreasing 
number of cis-isomers in the film the curvature of the 
film becomes smaller and the film eventually returns to 
its initial shape. This is happening in the experiments 
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FIG. 6: The same as in Fig. [S]for Tp — tv/A. 
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FIG. 7: The shape of the rectangular film with Aa > Ai 
illuminated by a linearly polarized light for ip — 0. 



of Refs. Q when the film is irradiated with light of 
the wavelength > 540 nm when the cis-isomeric state is 
depopulated. 



VI. CONCLUSION 

In conclusion, we presented a nonlinear elastic theory 
which describes an anisotropic bending of molecular films 
under the absorption of polarized light. Solving equa- 
tions of equilibrium, we showed that a change in the po- 
larization direction of light causes a corresponding change 
of the shape of the film. We also showed that the curva- 
ture of an irradiated film is a non-monotonic function of 
the extinction coefficient. 
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where 



dp \ dp dp 
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= 



Id d 
p dp dp 



(59) 



is the radial part of the Laplace operator. Inserting 
Egs. ipS)) into Eqs. (|^ - (P5)) . we obtain the boundary con- 
ditions at /9 = 1 for the zeroth harmonics in the form 



FIG. 8: The normalized mean curvature H/Hm as function 
of scaled extinction length ^//i. 
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APPENDIX 



^ = 0, 
dp 

and for the first Fourier harmonics in the form 



rfA " P' 



Ci+4 



1 - a 



Ci-p 



dCi 

dp 



^-aa = 
P 



4 \ 1 - o- , 

72 Ci + — ^ (4Ci- 



(62) 



(63) 



P^)-«. = 0,(64) 



We will look for solutions of Eqs. (PT|l - ([15|) in terms of 
Fourier series Introducing Eas. (P5|l into Eas. (|lT|) . 

(|42p and taking into account only the first two terms in 
the expansions (|48|l. we get 



A^o 



1 _d_ fdro dXo 
p d p \ d p dp 



1 ci Ci dTi 

2 dp dp 
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dp \p 



)} 



-/ 



dp 
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dp \p 



2pdp(\dp 



(56) 



(57) 



ri = , ^ - 



dri 
dp 



(65) 



As it is seen from Eqs. (j63| . ([64|) the amplitude of the 
first harmonics is proportional to the light-film interac- 
tion parameter Oq. Assuming that Oa < 1, we neglect all 
terms proportional a" with n > 2. Under this assump- 
tion one can neglect the last two terms in Eqs. ([5S|) . ([57|) 
and obtain 



^ pdpldpdpi 



A^ro 



1 1 d 

p dp 

Integrating each equation, we get 



(66) 



p2 



Ci 



^_d_ /drodCA _ 4 A dVo 
dp \ dp dp J p \ dp^ 



d / d(o dri 
dp \ dp dp 



dp^ 



= 0, (58) 



did dCo 
d p p d p dp 



drp dCo 
dp dp 



did drn , 9, 1 

d p p d p d p 2 



in/ 



dp 



1), 



=0, (68) 
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where the condition of regularity at the center d(^o/dp = 
and the boundary conditions ([50]) . ([5^ were used. 

Equations (|68p are simphfied by introducing the new 
variables 



1 dCo 11 dTQ 2 

pup A p dp 



(69) 



Then Eqs.dSHl) become 



■f 1^ 



— i^g)-ga=--{l--), 
d^ , , 3(1 -(t2) 2 „ 



(70) 



The boundary conditions at z = 1 become 



d^ 

2z^ + (l + cr)5-ai = 
a z 

a = . 



(71) 



Assuming that e = 3(1 — o'^)/8 < 1 is a small parameter 
we expand the functions g(z) and a{z) into scries 



9 = 90 + f^gi -\ — , 
a = OfQ + e ai + ■ • • . 

Inserting ^ into Eqs.dTD]) we get 

I 1 

^° • -mi^-go) - 9oaa = --{I- -) 
dz^ 8 z 



(72) 



^(zao) = 



■j^ (zgi) - .90 ai - 51 Qfo = 
j^izai)+ g^ = . 



(73) 



(74) 



In the same way the boundary conditions (|7ip aX z ~ 1 
can be expressed as 



^(^5o) = 



dan , 
2z-p+ 1 + a 5o-a» = 0, 
a 2; 

ao = , 



: ^(251) = 0, 

2z^ + (l + a)5i=0, 
ai = . 

From Eqs. (|73| and (|75|) one finds 

fl j f ( 1 - (7 z 

90 



, , , H In z 

l + (7 8 t2(l + cr) 2 



tto = 



(75) 



(76) 



(77) 
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